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Short gamma-ray bursts from SGR giant flares and 
neutron star mergers: two populations are better than one 
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ABSTRACT 

There is increasing evidence of a local population of short duration Gamma-ray Bursts 
(sGRB), but it remains to be seen whether this is a separate population to higher 
redshift bursts. Here we choose plausible Luminosity Functions (LF) for both neutron 
star binary mergers and giant flares from Soft Gamma Repeaters (SGR), and combined 
with theoretical and observed Galactic intrinsic rates we examine whether a single 
progenitor model can reproduce both the overall BATSE sGRB number counts and a 
local population, or whether a dual progenitor population is required. Though there are 
large uncertainties in the intrinsic rates, we find that at least a bimodal LF consisting 
of lower and higher luminosity populations is required to reproduce both the overall 
BATSE sGRB number counts and a local burst distribution. Furthermore, the best fit 
parameters of the lower luminosity population agree well with the known properties 
of SGR giant flares, and the predicted numbers are sufficient to account for previous 
estimates of the local sGRB population. 
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1 INTRODUCTION 

Results from the Burst and Transient Source Experiment 
(BATSE) onboard the Compton Gamma-Ray Observatory 
showed that Gamma-ray Bursts (GRBs) divide observa- 
tionally into two classes based primarily on their du- 
ration |Kouveliotou et all 1 19931 ): long GRBs have dura- 
tions > 2 seconds, and short GRBs sC 2 seconds. Short 
GRBs (sGRBs) seem to be associated with a variety 
of host ga laxies with no apparent restriction on galactic 
prop erties l|Prochaska et ahlboOrj ; iBerger 20(53; iLevan et al.l 
2008) , although host identification is not always triv- 
ial jLevan et al.l I2007T ). Additionally, a handful of recently 
detected sGRBs h ave localisations consist e nt with origins 
in nearby galaxies (Ipfek et alj 20061. 2007a |_ Frederiks et alj 



120071 ; iMazets et ail 120071 ; ILevan et al.1 120081). Over all, the 
Swift redshift distribution of sGR Bs (iBergerl 12007 ) peak s 
closer than that of long GRBs (jjakobsson et aL [2006), 
though there is evidence that some sGRBs may occur at 
higher redshifts (|Levan et al.l l2006al ). and that there may 
be a local popul a tion of underlum i nous l ong GRBs (e.g. 
Pian et al.l I2OO6I ; ISoderberg et al.1 120061; ILiang fc Zhand 



20061 ; iLiang et al.ll2007l : IChapman et al.ll2007h " 



The leading progenitor model for sGRBs is the merger 
of two compact objects, neutron star-neutron star (NS-NS) 
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or neutron star-black hole (|Nakarll2007f ) binaries. The Lu- 
minosity Function (LF) of BATSE sGRBs has been in- 
vestigated preyjously_^S i umi^a_sm 

tion ( e.g. ISchmidul200ll: lAndoll2004 ICuetta fc Piraiill2005l . 
2006; iHopman et al.1 120061 ) in order to determine the in- 
trinsic rate and most likely LF parameters. In a refine- 
ment to their previous work, iGuetta fc Piranl l|2006l ) noted 
that a second population of bursts may be necessary to 
explain some features of their model fits and the com- 
parison with Swift b ursts, particularly at lower redshifts. 
IHopman et al.l (|2006l ) considered both primordial and dy- 
namically formed NS binaries, and suggested that the early 
observed redshift distribution of s GRBs favoured dynami - 
cal formation. Further to that work. [Salvaterra et al.l (|2007l ) 
suggested that the more recent Swift cumulative redshift 
distribution is better encompassed by including both for- 
mation routes with different abundances above and below 
2 ~ 0.3. Recently, in an analysis of a large number of models 
of compact object merger scenarios from population s ynthe - 
sis models, lO'Shaughnessv. Belczvnski fc Kalogeral (|2008l ) 
have shown that the observed sGRB redshift distribution 
could be reproduced by a reasonable fraction of those mod- 
els, though this analysis was insensitive to the low end of 
the redshift distribution on w hich our work here is focused. 
iNakar. Gal- Yam fc F ox (2006) find the high rate of observed 
sGRBs within 1 Gpc to imply that a single population of 
NS binaries responsible for all sGRBs must be dominated 
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by long merger times, inconsistent with the observed NS 
binary population. However, they also point out that a non- 
unimodal luminosity function, such as produced by two sep- 
arate populations of progenitor, cannot be ruled out for 
sGRBs. 

There are indeed other possible progenitors for sGRBs. 
At much closer distances still, the initial spike in a giant 
flare from a Soft Gamma Repeater (SGR) in a relatively 
nearby galaxy would also appear as a sGRB. For example, 
the December 27th 2004 event from SGR 1806-20 would have 
been visible by BATSE out to « 50 Mpc jHurlev et al 1120051 ; 
IPalmer et all 120051 : iTavlor fc Granotl 120061 : iNakari l2007h . 
Thus it is entirely plausible that some fraction of sGRBs 
are extragalactic SGR giant flares. Several studies have es- 
timated theJikel^_coMri^utions of SGR flares to BATSE 
sGRBs. IPopov fc Sternl ()2006l ) estimate a rate of a few per- 
cent based on a lack of definite sGRB detections among 
the shortest BATSE GRBs consistent with locations within 
the Virgo cluster. Searches for hosts plaus i bly c onnected 
with six well-localised sGRBs (|Nakar et alj 120061 ) suggest 
a rate of less than 15%, and in a comparison of the spec- 
tra of a limited sample of the bright est BATSE sGRBs, 
Lazzati. Ghirlanda fc Ghisellinil |2005) conclude only 4%. 



Palmer^alJ M), M ^ta^l lack of events from 



the Virgo cluster, find a rate of less than 5%, though point 
out that the LF of SGR giant flares may extend to much 
larger luminosities, such as sugg ested bv lEichlerl (|2002D . IOfekl 
l|2007h points out that the fraction cannot be less than ~ 1% 
without being inconsistent with the observed Galactic SGR 
giant flare rate, and calculate an upper limit of 16% (95% 
confidence limits) based on a conservative measure of proba- 
ble IPN sGRB coincidences with bright star forming galaxies 
within 20 Mpc. This limit is sensitive to their estimate of 
the completion of the galaxy sample and may be higher still. 

Previously, using the full sample of BATSE sGRBs lo- 
calised to better than 10 degrees, we demonstrated that be- 
tween 6 and 12 per cent of BATSE sGR Bs were correlated o n 
the sky with galaxies within ~ 28 Mpc l|Tanvir et al.| [2005). 
and we have now extended this work out to ~ 155 Mpc. 
Our analysis was based purely on burst /galaxy distribution 
correlations and unbiased with regards to burst brightness 
or other assumptions, though our later work showed that 
this correlation is dependent mostly on large scale struc- 
ture on the sky rather tha n individual burst/galaxy pair- 
ings l|Chapman et al.l 120071 ). The main question we address 
in this paper is whether a nearby population (z ^ 0.03) of 
this magnitude may be produced by a suitable LF describing 
a single progenitor population, or whether it is necessary to 
include an intrinsically lower luminosity population as well. 

Here we attempt to answer this question by consider- 
ing first single, and then dual population LFs. The intrinsic 
rates in the models will be assumed from both the observed 
Galactic SGR flare rates and the modelled NS-NS merger 
rates in order to investigate the LF parameters. Obviously 
there are significant uncertainties in these rates: the Galac- 
tic giant flare rate in particular is estimated from only 3 
observed events. Regardless of these uncertainties and the 
exact form of luminosity functions chosen, we find that a sin- 
gle progenitor population described by a unimodal (i.e. with 
a single peak or knee) LF cannot produce sufficient local 



events, whereas a dual population reproduces the likely lo- 
cal sGRB distribution as well as the overall number countfl 



2 METHODS 

The number of sGRBs, N, observed above a threshold p in 
time T and solid angle Q is given by Equation[T] where <E>(L) 
is the sGRB LF, Rgrb(z) is the comoving event rate density 
at redshift z, dV(z)/dz is the comoving volume element at 
z and z max for a burst of luminosity L is determined by the 
detector flux threshold and the luminosity distance of the 
event. 



N(>p) 



4% 



Rgrb(z) dV(z) 



m Jo 1 + z dz 

(1) 

We are of course dealing with detector l i mited and 
not bolometric luminosities. Following ISchmidtJ (|200lh and 
iGuetta fc Piranl (|2005h we assume a constant median spec- 
tral index of -1.1 in the BATSE energy range of 50-300 keV 
to derive a simplified K correction and conversion to photon 
flux. 



2.1 Intrinsic rates 

The sGRB rate per unit volume, -Rgrb(z) is given by Equa- 
tion [5] where Nqkb is the number of sGRBs per progenitor, 
Pprogenitor is the intrinsic (z = 0) progenitor formation rate 
and F(z) describes the volume evolution of this rate with z. 



RGRB{z) = N G RB X /Oprogcnitor X F(z) MpC 



(2) 



For NS-NS mergers, a burst is produced only once at 
merger, and we therefore assume ATqrb = 1. This is of course 
an upper limit: any beaming of sGRBs, or a GRB produc- 
tion efficiency per merger of less than 100%, would effec- 
tively reduce this number, and reduce the number of bursts 
observable from the NS-NS merger population. This limit 
is therefore conservative in the sense that it maximises the 
possible fraction of bursts produced by mergers in our anal- 
ysis. The intrinsic NS-NS merger rate is taken as 10 -5 yr _1 
per Milky Way equival ent galaxy ( Star Formation Rate, 
SFR « 4M Q yr" 1 , e.g. iDiehl et all d2006l)) from the pop- 
ulation synthesis models of iKalogera et al (|2007T ). Mergers, 
of course, occur some time after the formation of the bi- 
nary itself. Thus the merger rate at redshift z, is depen- 
dent not on the SFR at the same z, but on the earlier SFR 
at higher redshift. F(z) is therefore given by the convolu- 
tion of the SFR as a function of redshift with a distribu- 
tion of delay times fro m binary formatio n to m erger. The 
population syntheses of iBelczvnski et al] (|2006h suggest a 
merger delay time (formation plus coalescence) distribution 
dP m /d(log(t)) ~ constant (= dP m /dt oc 1/t) between 10 7 
and 10 10 years, with a narrow peak at the very lowest times, 
and we thus assume a delay time probability distribution 
where dP m /d(log(t)) is flat between 10 7 and 10 10 years and 
zero outside this range, for simplicity and comparison with 



1 Throughout this paper we assume a flat cosmology with Ho 
71km s _1 Mpc -1 , n M = 0.27 and U A = 0.73 
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previous LF analyses. We note, however, that using a de- 
lay model including a narrow early 'spike' (with an order 
of magnitude higher value between 15 and 30 Myr) makes 
little difference to the derived LF parameters as can be seen 
from some examples in Table [1] 

SFR as a fu nction of z is parameteri sed according to 
the SF2 model of lPorciani fc Madaul d200lf). normalised to a 
local SFR of 1.3 x 10" 2 M^vr~'Mpc 3 l|Gallego et al.lll995l ) 
as given in Equation [3] 

SFR(z) = 1.3 x 1(T 2 ( ^l 6 ^ ) M ° yr_1 MpC ~ 3 (3) 

An alternative analysis is that the merger rate should 
be proportional to Stellar Mass Density (SMD), which must 
be representative of star formation history. We therefore also 
investigate merger rates which follow a simple single expo- 
nential fit to the SMD out to z ~ 5 derived from the FORS 
deep field jProrv et al.ll2005T l as: 



and 



10 47 erg s 1 |Tanaka et al.ll2007f ) (including a correc- 



SMD(z) = 10 8 ' 75 exp(- 



ln(2)z) M Mpc" 



(4) 



Over the last 30 years of observations, there have been 
3 giant flares from 4 known SGRs in the Milky Way and 
Magellanic Clouds. The observed local rate of giant flares per 
Galactic SGR is therefore w 3 x 10~ 2 yr -1 , and their short 
active lifetimes of ~ 10 years l|Duncan fc Thompsonlll992l ; 
iKouveliotou et al.|[l99Sl ; lKouveliotoulll999l ) imply A^rm ~ 
300 in the SGR case. Magnetars are commonly believed to 
form in a fraction of core collapse supernovae, and hence 
their formation should follow the SFR as a function of z. 
Given the the association of the 4 known SGRs with young 
stellar populations, this therefore implies a formation rate 
via core collapse supernovae of 4 x 10~ 4 yr _1 . 

However, it is also plausible that magnetars may form 
via the Accretion Induced Collapse (AIC) of White Dwarf 
(WD) binaries which contain at least one sufficie ntly mas- 
sive and magnetized member (|Levan et al . 2006b). In older 
galaxies with relatively little star formation, this would be 
the dominant formation route and therefore makes it pos- 
sible for SGRs to be associated with all types of galax- 
ies, not just those with a relatively high SFR. Following 
iLevan et al.l (|2006bL the rate of magnetar formation via 
WD- WD mergers in a Milky Way equivalent galaxy is es- 
timated as 3 x 10 -4 yr _1 . We therefore assume F(z) for 
SGRs follows both SFR(z) for magnetar production from 
supernovae and either the delayed SFR or SMD to allow for 
production by WD binary mergers. 

2.2 Luminosity functions 

Luminosity functions for SGR giant flares and NS-NS merg- 
ers are not well known. A lognormal LF approximates the 
shap e of the theoretical NS-NS merg er luminosity distribu- 
tion (|Rosswog fc Ramirez- Ruizl 12003), but other functional 
forms may be equally valid: for example iGuetta fc Piranl 
(2005) assumed a broken power law for their LF calculations, 
and the luminosities of many other astrono mical population s 
are well described by a Schechter function (|Schechterlll97r3) . 

Given only 3 events, it is not possible to constrain 
the SGR giant flare LF to any great degree. The 3 ob- 
served events have peak luminosities of ~ 10 44 , ~ 10 46 



tio n for the lower dis tance estimate of SGR 1806-20 found 
bv lBibbv etaLI (|2008h V The more common short duration 
bursts from SGRs, with luminosities up to 10 41 erg s _1 seem 
to follow a power law distribution in energy, dN oc E r dE 
where 7 ~ 1.4 - 1.8 (|Cheng et al.lll996l ; IGogiis et alJl200Ch 
similar to that found in earthquakes and solar flares. Inter- 
mediate bursts with energies and luminosities between the 
short bursts and giant flares are also seen, and it is possible 
therefore that this distribution continues to higher energies 
and includes the g iant flares themselves, particularly since 
IGogiis et al.l (|2000T l found no evidence for a hig h energy cut- 
off in their work. However, ICheng et all l|l996h did find evi- 
dence of a cutoff around 5 x 10 41 erg, and furthermore the 
intermediate bursts are generally seen following giant flares 
and may be some form of aftershock rather than represent- 
ing part of a continuous spectrum of flare activity. Theory 
suggests that the common bursts are produced by the re- 
lease of magnetic energy gated by a small scale fracturing of 
the crust sufficient only to relieve crustal stresses, whereas 
the giant flares are the result of large scale cracking suffi- 
cient to allow external field reconfiguration to a new equi- 
librium state ^Thompson fc Duncan|[l993l . |l995T ). Assuming 
the latter is a physically distinct process discontinuous (in 
terms of energy release) from the short bursts, then it must 
have some minimum energy release, and a maximum de- 
fined by the total destruction of the external field via the 
Flowers- Ruderman instability (jFlowers fe Rudermanl [19771 ) 
where entire h emispheres of the magnetar flip with respect 
to each other (Eichler 2002). Having only the 3 observed 
events to go on, a lognormal LF is once again plausible for 
giant flare luminosities. The possibility of a continuous lu- 
minosity distribution between the short, intermediate and 
giant flares is not ruled out however, and we therefore also 
consider a single power law LF as well. 

To summarise, we consider the possibility that short 
GRBs may be produced via two different progenitor routes, 
both NS-NS mergers and SGR giant flares, each population 
with intrinsically different luminosities. The forms chosen 
for the luminosity functions examined are as follows: 
1. Lognormal distribution 



dlog L 
2. Schechter function 



dN /-(logL-logio) 
oc exp ( 



2<r 2 



(5) 



^ oc (yj^j exp(-L/L ) , L > L min (6) 



3. Power Law 



dN 
~dL 



<L<L 



(7) 



where L min — 10 42 erg s _1 for normalisation and con- 
vergence of the Schechter function (see Appendix A for dis- 
cussion of the limited effect of the choice of Lmin)- The 
Power Law distribution is normalised to the observed Galac- 
tic rate between LaFmin = 10 44 erg s _1 and Lo, but the 
distribution is analysed down to L m i n to investigate the pos- 
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sible extension of the power law to lower luminosity flares. 
Lo, and a or a are the free parameters to be estimated. 

2.3 Constraining the models 

The Cmax/Cmin table fro m the current BATSE cata- 
logue l|Paciesas et alj Il999l ) provides peak count rate for 
bursts in units of the threshold count rate. Not all bursts 
are included and in addition the BATSE threshold was var- 
ied historically. Therefore in order to analyse a consistent 
set of bursts we restricted the table to only those sGRBs 
recorded when the 64ms timescale threshold was set to 5.5<r 
above background in at least 2 detectors in the 50 — 300 keV 
range. The all sky equivalent period (including correction 
for BATSE's sky coverage) this represents is estimated as 
~ 1.8 years. 

We then examined the differential distributions of pre- 
dicted overall counts first from various single, and then com- 
bined populations of burst progenitor. By varying the pa- 
rameters of the chosen luminosity functions, we compared 
the predicted overall counts (dN/dp) to the observed differ- 
ential distribution from the Cmax/Cmin table. For each set 
of LF parameters, the redshift distribution of sGRBs was 
calculated, and the nearby distribution s compared with th e 
observed correlated distributions from lTanvir et all (jj2005) . 
Note that we use an extended version of our previous cor- 
relation analysis out to 155 Mpc, and use the correlations 
measured against galaxies in concen tric shells (as opposed 
to spheres) of recession velocity (see IChapman et all 12007! ) 
in order to obtain a local differential distribution for the 
model fitting, x 2 minimisation was then used to optimise 
the LF parameters by fitting simultaneously to the overall 
count rate and the local distributions. We assumed a Pois- 
sonian error distribution on the overall count rate, whereas 
we used the explicit Monte Carlo derived error distribution 
on the local correlated fraction (the error distributions from 
the Monte Carlo simulations closely follow a normal dis- 
tribution, even at low correla tion levels since the function 
defined in iTanvir et alj l|2005h is equally sensitive to anti- 
correlation giving rise to negative percentage correlations 
in those situations). Note that the greater number of data 
points in the number count fits means that the combined 
\ 2 values are dominated by the goodness of fit to the count 
rate distribution. To explicitly ask whether any of our chosen 
single LFs can remain consistent with the BATSE number 
counts while being forced to produce a local distribution of 
bursts, we also find the best fits constrained by the corre- 
lated fraction alone. 

In order to check the plausibility and consistency of 
the best fit models, we further compared the derived red- 
shift distribution to that of sGRBs observed by Swift. We 
caution that this sample is neither uniformly selected nor 
complete. Previous studies have an alysed the early Swift 



distributions (e.g. Guctta fc Piranl [2006; Hopman et all 



20061; iNakar. Gal- Yam fc Foxl l2006t ISalvaterra et all 120071 : 



O'Shaughnessv. Belczvnski fc Kalogeral 20081 ). and it is 
clearly useful to compare our models to the current best 
known redshift distribution in order to check that the predic- 
tions are not unrealistic. We stress that the Swift distribu- 
tion was not part of the statistical analysis. sGRB redshifts 
have so far only been found from host galaxy associations, 
the identification of which is not always unambiguous. Fur- 



Table 1. Results of single population Luminosity Functions, 
presented in order of decreasing overall goodness of fit (i.e. in- 
creasing overall \ 2 /dof). The LFs follow merger delay time (for- 
mation plus coalescence) distributions either flat in log space 
(dP m I d(log(t)) = constant), or with a narrow spike at early 
times, or the SMD profile of Equation [4] The number of degrees 
of freedom (dof) for the Cmax/Cmin and overall distributions 
are 22 and 26 respectively. Iq is in units of log(erg s — x ), a in dex 
and a is dimensionlcss. 



NS Merger 


Parameters 


Local 


Cmax /Cmin 


Overall 


LF 


(lo = 


Ogio) 


x 2 


X 2 /dof 


X 2 /dof 


Schechter 


k = 


51.75 


11.93 


1.01 


1.31 


(flat) 


a = 


1.25 








Schechter 


lo = 


51.8 


11.93 


1.03 


1.33 


(spike) 


a = 


1.25 








Schechter 


lo = 


50.45 


11.87 


1.05 


1.34 


(SMD) 


a = 


= 0.9 








Lognorm 


lo = 


48.9 


11.82 


1.09 


1.38 


(SMD) 


a = 


0.75 








Lognorm 


lo = 


46.4 


11.88 


1.18 


1.46 


(flat) 


a — 


1.5 








Lognorm 


lo = 


46.6 


11.89 


1.19 


1.47 


(spike) 


a = 


1.45 









thermore, even the classification of some bursts as either 
short or long is controversial since their durations change 
substantially depending on whether or not emission from 
the long-soft tails (seen in a number of bursts) is included. 
Nevertheless, about a dozen probable short-hard bursts 
have reasonably secure redshifts. Specifically we include 
the following 10 sGRBs: GRBs 050509B, 05 0724, 051221a , 
060801, 061006, 061201, 061210, 061217 fseelBergerl (I2007T ) 



and refe rences therein) 07071 4B (jGraham et alll2007l) and 



071227 (|D'Avanzo et al.| [2007). In order to produce the pre- 
dicted Swift redshift distribution, the Swift BAT thr eshold 
for sGRBs was assumed to be twice that of BATSE l|Bandl 
I2006D . 



3 RESULTS 

Table [1] lists the best fit parameters found from fitting dis- 
tributions produced by single population NS merger LFs 
simultaneously to both the overall number counts and the 
local population as described above. The table is ordered 
in decreasing overall goodness of fit (i.e. increasing com- 
bined x 2 /dof). As mentioned previously, the combined \ 2 
is dominated by the fit to the overall BATSE number counts 
and, as expected, all our chosen single population LFs pro- 
duce good fits to the Cmax/Cmin data leading to acceptable 
overall fits as measured by the combined \ 2 ■ However, none 
of the single progenitor population LFs reproduce the local 
burst population expected from the correlation results: for 
example Figure [T] shows the results from a single Schechter 
function LF which can be seen to produce effectively no 
sGRBs within 300 Mpc. To enable a quantitative compari- 
son with later results, column 3 in TableQ]lists the \ 2 results 
considering the fit to the four data points of the local dis- 
tribution alone. Note that since the fit is constrained by the 
overall data set, the local x 2 value is quoted unreduced in 
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Total no. of bursts 



□ No. of bursts from Tonvir et ol, 2005 



Table 2. Results of single population Luminosity Functions con- 
strained to fit the local distribution, presented in order of decreas- 
ing goodness of fit. Details as for Table [T] 



200 300 
Distance (Mpc) 



i» 0.8 




Total no. of bursts 
Swift observed bursts 



Figure 1. Burst distributions from the best fit merger single 
population Schechter function LF following a dP m / d(\og(t)) = 
constant merger time delay distribution. Top panel shows pre- 
dicted sGRB distribution wit hin 500 Mpc comp ared to the lo- 
cal burst fraction measured in lTanvir et al.l l|2005h . bottom panel 
shows the predicted burst distribution out to z = 3 normalised 
and compared to the Swift distribution discussed in the text. 



this and all subsequent Tables since the precise number of 
degrees of freedom relevant to this subset alone is difficult 
to estimate formally. 

In order to ascertain whether a single merger population 
can produce the local bursts and remain consistent with the 
Cmax/Cmin data, we then fit single LF populations to the 
local distribution alone, with no constraints placed on good- 
ness of fit to the overall number counts. As can be seen from 
Table [5J single Schechter function LFs can produce a local 
population, but the associated number count distribution is 
an extremely poor match to the C max / C m in data. 

Of course, these results represent the best possible re- 
production of the local bursts (in terms of minimum \ 2 val- 
ues), and it maybe that there exist poorer fits to the local 
population which are nevertheless better fits to the overall 
number counts. Figure[2]shows the \ 2 contours (individually 
for fits to both the overall number counts and the local dis- 
tribution) for the single population Schechter function LF 
(dP m / d(\og(t)) — constant merger delay time distribution) 
from Table Q] From the Figure it can be seen that accept- 
able fits to the overall number count distribution (to a sig- 
nificance level of 0.999) are found only in a narrow band of 
the LF parameter space, well separated from even the 0.99 
significance level of the local distribution fits (which repre- 
sents effectively no local distribution given the size of the 
errors on the correlation points shown in the top panel of 



NS Merger 


Parameters 


Local 


Cmax 1 Cjnin 


Overall 


LF 


(lo=i 


ogL ) 


x 2 


X 2 /dof 


X 2 /dof 


Schechter 


lo = 


51.1 


0.76 


> 100 


> 100 


(flat) 


a = 


2.2 








Schechter 


lo = 


51.0 


0.76 


> 100 


> 100 


(spike) 


a = 


2.2 








Schechter 


lo = 


53.0 


1.08 


> 100 


> 100 


(SMD) 


a = 


2.15 








Lognorm 


lo = 


43.1 


5.30 


> 100 


> 100 


(spike) 


cr = 


1.35 








Lognorm 


lo = 


43.1 


5.38 


> 100 


> 100 


(flat) 


(7 = 


1.3 








Lognorm 


lo = 


43.0 


6.43 


> 100 


> 100 


(SMD) 


a = 


1.4 










Figure 2. \ 2 contours for the single population Schechter func- 
tion LF (dP m /d(log(t)) = constant merger delay time distribu- 
tion) from Table [T] The solid lines are the contours for the fits 
to the overall number counts, and the dotted lines are contours 
for the fits to the local distribution. Contours are plotted at 0.68, 
0.9, 0.99 and 0.999 significance levels. For the sake of clarity, only 
the outermost 0.999 contour is labelled for the number count fit. 



Figure [T] for example). Hence the possibility of this single 
Schechter function LF to reproduce the local distribution 
while remaining consistent with the overall number counts 
can be rejected with greater than 99.9% confidence. Similar 
results with equivalent levels of rejection are found for the 
other single LF models. 

The inability of the local constraint to produce a distri- 
bution which fits the number counts is effectively a conse- 
quence of the intrinsic sGRB rate calculated in Equation [2] 
from the assumed merger rates: not enough bursts in total 
can be produced. The fit to the Cmax/Cmin data can be 
improved to reasonable \ 2 levels by increasing the intrinsic 
merger rate by a large factor (^ 500), but the overall red- 
shift distribution produced as a consequence is extremely 
unrealistic, with all bursts produced within z ~ 0.1. 

In contrast, Table [3] shows best fit LF parameters for 
various combinations of dual NS merger and SGR giant flare 
luminosity function models, along with their respective min- 
imum x 2 values. As can be seen from Table [2] and the ex- 
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Total no. of bursts 

Bursts due to SGR giant flares - 

Bursts due to NS-NS mergers ; 

No. of bursts from Tanvir et al, 2005- 



200 300 
Distance (Mpc) 



0.3 



0.2 




0.0 



Total no. of bursts 

Bursts due to SGR giant flares - 

Bursts due to NS-NS mergers ; 

□ No. of bursts from Tanvir est -air 20"05 : 



200 300 
Distance (Mpc) 




Total no. of bursts 
Bursts due to SGR giant flares- 
Bursts due to NS-NS mergers 
Swift observed bursts 
Swift observed (inc. 050906) 



3 0. 




Total no. of bursts 
Bursts due to SGR giant flares- 
Bursts due to NS-NS mergers 
Swift observed bursts 
Swift observed (inc. 050906) 



Figure 3. Burst distributions from dual lognormal LF (follow- 
ing dP m / d(\og(t)) = constant merger time delay distribution) 
populations. Panel Details as for Figure [l] 



amples of Figures [3] and [H the local bursts are only ever 
reproduced by the lower luminosity LF. By minimising the 
combined \ 2 values, all the dual LFs tested reproduced the 
local distribution well while retaining overall number count 
fits comparable to those of the single LFs. Furthermore, the 
best fit LF parameters of the dual models are reasonable, 
and the overall redshift distribution is much more realistic. 

For example, a dual lognormal LF, with merger rates 
following either a delayed merger model (Figure [3} or the 
SMD model of Equation [4] (Figure [4}, produces a good fit 
to the expected local population while remaining consistent 
with the early Swift redshift distribution. The upper panels 
of Figures [3] and [4] show the comparison of these models 
to the local sGRB distribution determined by our BATSE 
cross-correlation analysis, and are typical in that all the dual 
populations reproduce this local population well. Since these 
data were used to constrain the fit, a good agreement is to 
be expected, but it is still interesting to note that the merger 
population contributes only a small fraction to these local 
bursts. The lower panels show the overall predicted redshift 
distribution. 

As mentioned before, the intrinsic Galactic rates used to 
normalise the LFs are not well constrained. Hence in Table[3] 
we also show the results of varying the intrinsic SGR flare 
rate up and down by an order of magnitude for the dual log- 
norm (SMD) fit of Figured The production of a local sGRB 
population is robust against this change, and the overall fit 
remains good. As may be expected, an increase in the in- 
trinsic flare rate leads to the best fit SGR LF being moved 



Figure 4. Burst distributions from dual lognormal LF (following 
SMD) populations. Panel Details as for Figure [T] 

down in luminosity, thus removing a greater fraction of the 
total flares from observability. Likewise, a lower intrinsic rate 
generates a higher (and narrower) LF distribution, though 
in both cases the LF parameters remain entirely plausible. 

Figure [5] shows the best fit LFs and associated contours 
of x' 2 with respect to Ln for the dual population from Fig- 
ure [4] Despite the uncertainties in the underlying Galactic 
rates of the models, the best fit parameters obtained for 
this and the other dual LFs are plausible given the known 
properties of SGR giant flares and classic sGRB luminosi- 
ties. We note that the slopes of the SGR flare power law 
LFs obtained (1.25 - 1.35) are shallower than the slopes 
foun d for ordinary SGR burst fluence di stributions (1.4 - 
1.8) (|Cheng et al.lll996l : iGogiis et alJl200Ch . 



4 DISCUSSION 

The lower panels of Figures [3] and 3] imply that Swift should 
have triggered on about one SGR flare to date (this would 
rise by a factor of ~ 2 if the redshift completeness for such 
flares were greater than for sGRBs as a whole, as is likely 
given that low-redshift host galaxies are easily identified). 
We note that a possible candidate is GRB 05 0906, which 
may have originated in a galaxy at « 130 Mpc (|Levan et al.l 
2008), and the preliminary Swift redshift distributions in 
Figures and [3] are plotted both including and excluding 
this burst. 

There are two further recent sGRB events which 
are candidate extragalactic SGR flares, though neither 
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Table 3. Results of dual population Luminosity Functions, presented in order of decreasing goodness of fit (i.e. increasing overall \ 2 /dof). 
The LFs follow merger delay time (formation plus coalescence) distributions either flat in log space (dP m / d(log(t)) = constant) or the 
SMD profile of Equation|4] Also shown are two results normalised using order of magnitude different observed Galactic (MW) rates. The 
number of degrees of freedom (dof) for the C m ax/C m i n and overall distributions are 20 and 24 respectively. lo is in units of log(erg s — 1 ), 
cr in dex and a is dimensionless. 



NS Merger LF Parameters SGR Giant Flare LF Parameters Local C ma x/C m i n Overall 
LF (jo = logic,) LF (Jos log Lo) X 2 X 2 /dof x 2 /dof 



Schechter 


lo 


= 52.3 


Power law 


(flat) 


a 


= 1.3 


(flat) 


Schechter 


h 


= 52.3 


Lognorm 


(flat) 


a 


= 1.3 


(flat) 


Lognorm 


lo-- 


= 48.35 


Lognorm 


(SMD) 


a 


= 1.05 


(SMD) 


Lognorm 


lo 


= 47.2 


Power law 


(flat) 


<T 


= 1.2 


(flat) 


Lognorm 


lo-- 


= 47.05 


Lognorm 


(flat) 


a 


= 1.2 


(flat) 



lo = 46.7 


2.03 


1.15 


1.04 


a = 1.25 








l = 45.2 


1.45 


1.20 


1.06 


cr = 0.6 








lo = 45.3 


1.66 


1.31 


1.16 


a = 0.55 








l = 46.7 


2.06 


1.69 


1.49 


a = 1.35 








Z = 45.2 


1.55 


1.72 


1.50 


cr = 0.6 









Lognorm l = 48.6 Lognorm l = 44.1 1.57 1.36 1.20 

(SMD) a = 0.9 (SMD) cr = 0.8 

(10 x MW) 

Lognorm l = 48.6 Lognorm Z = 46.3 3.13 1.28 1.20 

(SMD) cr = 0.9 (SMD) a = 0.2 

(0.1 x MW) 



triggered Swift: GRB 051103 whose IP N error box in- 
clude s the outskirts of M81 at 3.5 Mpc |Golcnctski i et all 
120051 ). and GRB 070201 whose error box similarl y over- 
laps a spiral arm of M31 at only ^ 0.77 Mpc |Perlevl 
120071 ; iPal'Shinl 120071 : iMazets et aLll2007f). Both have char- 
acteristics o f SGR giant flares dFrederiks et al.l 120071 ; 
IMazets et ail 120071 ; lOfek et all l2007bl ). and furthermore 
the non-dete ction of gravitational waves by LI GO from 
GRB 070201 l|LIGO Scientific Collaboration! I2007T I excludes 
a merger progenitor within M31 with > 99% confidence. If 
both these events were due to extragalactic SGRs then this 
brings to three the number of giant flares with peak lumi- 
nosity > 10 47 erg s _1 seen in just a few years. 



iLevan et al.l (|2008l ) estimated that a Galactic SGR gi- 
ant flare rate of ~ 0.5 x 10 -4 yr _1 would be sufficient to 
produce ~ 10 extragalactic flares within a sphere of radius 
100 Mpc. Using a power law LF (constrained by a search 
for positional coincidences between galaxies within 20 Mpc 
and t he IPN error boxes of a sample of 47 sGRBs), lOfekl 
<|2007l ) estimated the rate of extragalactic flares with en- 
ergy > 3.7 x 10 46 erg (th e energy of the 2004 SGR1806-20 
event <|Hurlev et al.ll2005l )) to be ~ 0.5 x 10 -4 yr" 1 per SGR, 
and the 95% confidence lower limit of the Galactic rate to 
be 2 x 10~ 4 yr^ 1 per SGR. Our analysis estimates the rate 



1 x 10~ 4 yr" 1 per SGR. We es timate 



these two values at 
the SFR of galaxies within 5 Mpc listed b y Ofek ( 2007) (with 
revised distance estimates (|Karachentsev et al. |2004| )) to be 
about 22 x that of the Milky Way. Adopting our predicted 
(lognorm following SMD) flare rate, the probability of ob- 
serving two (one) or more such flares within this volume 
during the 17 years of IPN3 observation is 1% (14%). This 
indicates we have been witness to a rather rare coincidence, 



and is perhaps suggestive that not both GRB 051103 and 
GRB 070201 are SGR flares. 



5 SUMMARY AND CONCLUSIONS 

We have examined a selection of plausible Luminosity Func- 
tions, singly and in combination, for both neutron star merg- 
ers and SGR giant flares as progenitors of short Gamma-ray 
Bursts. Assuming observed and theoretical Galactic intrin- 
sic rates, merger delay time distributions, Star Formation 
Rate and Stellar Mass Density parameterisations, we ex- 
clude both lognormal and Schechter type LFs for a single 
NS merger population of progenitor as being unable to pro- 
duce a nearby sGRB population while remaining consistent 
with overall BATSE number counts. Indeed, given that even 
a Schechter function (dominated by low luminosity events) 
cannot reproduce the likely local population, it is hard to 
conceive of any unimodal LF which could and still be con- 
sistent with the higher redshift distribution. We suggest that 
at least a bimodal LF, and therefore likely a dual popula- 
tion model, is necessary to account for the local population. 
Given the uncertainties in the intrinsic rates assumed, we 
cannot sensibly choose between the LF combinations, but we 
point out that the best fit LF parameters in all dual popula- 
tions considered are in reasonable agreement with the known 
properties of SGR giant flares and classic sGRBs, even when 
the intrinsic rate of Galactic SGR flares is varied by an order 
of magnitude in either direction. To put this another way, as 
is well known a single population Luminosity Function pro- 
vides a good fit to overall BATSE number counts, but we 
find that a separate, lower luminosity population of progen- 
itors is both required, and is sufficient, to reproduce a local 
sGRB population. Furthermore, the properties of this pop- 
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Figure 5. Best fit dual population LFs from Figure [4] The LFs 
(top panel: dotted line SGR giant flares, dashed line mergers) 
are lognormal with intrinsic merger rate components following 
the SMD model of Equation [4] The bottom panel shows contours 
of x 2 m l°g(io) space. Contours shown represent 0.68, 0.9 and 
0.99 confidence limits with the minimum \ 2 value plotted as an 
asterisk. 

ulation are in agreement with those observed from Galactic 
SGR giant flares. 
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APPENDIX A: THE LOW LUMINOSITY 
CUTOFF, Lmin, IN THE SCHECHTER 
LUMINOSITY FUNCTIONS 

It may be thought that the precise choice of the lower 
luminosity cutoff, L m in, necessary for convergence of the 
Schechter type Luminosity Functions has a significant ef- 
fect on the LF parameters found, and further on the abil- 
ity of that LF to reproduce the local population while 
remaining consistent with the overall number counts. We 
therefore tested a range of low luminosity cutoffs (L m i n = 
10 40 erg s" 1 to 10 46 erg s" 1 , and found that the best fit 
LF parameters are relatively insensitive to the chosen Lmin 
as shown in Table IA1I Furthermore, the separation in LF 
parameter space is maintained between the best fits to the 
number counts and the best fits to the local distribution 
regardless of choice of L m i n as shown for example in Fig- 
ures [Si] and [Ml 



Table Al. Results of single population Schechter Luminosity 
Functions with varying values of L m i n . Details as for Table [T] of 
main text. 
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LF Parameters 
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l = 52.00 
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a = 1.36 
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Figure Al. \ 2 contours for the single population Schechter func- 
tion LF as in Figure[5]of the main text, with L m i„ = 10 40 erg s — 1 . 
Details as for Figure [5] of main text. 
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Figure A2. \ 2 contours for the single population Schechter func- 
tion LF as in Figurc[5]of the main text, with L m i n = 10 46 erg s _1 . 
Details as for Figure [5] of main text. 



